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The aim of the present paper is to propose a new type of Monte-Carlo 
method which enables us to study the excited state of many-body systems. 
qq ■ It consists of the method of continued fractions (MCF) and a Monte-Carlo 

random walk for multi-dimensional integration. The convergence of the MCF 
and its accuracy are studied in detail for a one-dimensional model. As an 
application of this method to a realistic problem, we study the ground state 
of 4 He by using the Volkov potential. It is shown that our results are in good 
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agreement with those of previous works obtained through different types of 
formalism. This implies that our new method works effectively in the study 
of many-body nuclear systems. 

I. INTRODUCTION 



The exact treatment of nuclear many-body systems is an important benchmark for inves- 
tigating the dynamics of multi-nucleon systems. It is a testing ground of potential models 
and equations of motion, such as those involving multi-nucleon interaction, relativistic dy- 
namics, and possible modification of nucleon properties inside nuclei. 

At present, various methods to solve many-body problems have been developed. One 
of the standard methods is the Faddeev or Faddeev-Yakubovsky (FY) equation. The FY 
equation has been applied to the investigation of the ground state of 4 He using a realistic 
nuclear potential. A common problem in attempts to extend the 'exact' study of many- 
body systems is the rapid increase of the number of channels and dimensions as the number 



of particles increases. One way to overcome this problem is to diagonalize the Hamiltonian 
in a subspace spanned by appropriate bases chosen using a stochastic method. The other 
way is the direct calculation of the matrix elements using a Monte-Carlo method. In the 
Monte-Carlo method, one can perform an effective mult i- dimensional integration by taking 
the importance sampling from the essential part of the phase space. 

Extensive works have been carried out by using a variational Monte-Carlo method (VMC) 
§ and a Green's function Monte-Carlo method (GFMC). |,§ VMC and GFMC have been 
applied to the study of the ground state properties of nuclei with mass number A < 7 us- 
ing realistic nuclear potential. || The bound state energies and also the electro-magnetic 
properties of few-nucleon systems have been studied extensively. In VMC, we can obtain 
the upper bound of the ground state energy for which the matrix elements are evaluated by 
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Monte-Carlo integration. Simultaneously, we obtain an analytic form of the trial function 
which describes the ground state quite well. In GFMC, we start the calculation with a trial 
function obtained by VMC, and the imaginary-time evolution e~ rH extracts the improved 
ground state wave function from the trial function. Here, the imaginary-time evolution is 
carried out by using the quantum Monte-Carlo random walk method. Due to the e~ rH pro- 
jection, we only obtain the properties for the lowest state of the system with a given quantum 
number. Excited states in the diffusion Monte-Carlo method can be investigated using a 
diagonalization procedure in the large space shell-model calculation with an auxiliary-field 
Monte-Carlo approach |j],[7| and in one- and two-dimensional examples of GFMC. || For 
the purpose of conducting systematic studies in few-nucleon systems, it is necessary to in- 
vestigate excited states as well as the ground state. Therefore we should develop a method 
for this purpose in the context of few-nucleon physics. 

In this paper, we propose a new Monte-Carlo method to study the bound state prob- 
lems in few-nucleon systems. Here, our purpose is to find an effective method to extract 
eigenfunctions without referring to the imaginary-time evolution method, but still using the 
Monte-Carlo random walk. Therefore, the method is not restricted to the ground state 
problem. We use the method of continued fractions (MCF) as a central formalism in our 
new calculation method. MCF was originally proposed by J. Horacek and T. Sasakawa 
as a fast calculational method to evaluate T-matrix elements and was applied to the 
three-nucleon system. The wave functions which appear in the MCF iteration are generated 
so as to be orthogonal to the lower order ones. Then the matrix elements calculated by 
those wave functions rapidly converge. Moreover, this method is simple and allows us to use 
any function as a starting trial wave function and to adopt a nuclear potential including a 
many-body potential. 

The direct application of MCF with Monte-Carlo integration, however, causes a few 
serious problems: First, the original form of the MCF iteration formula 0] is not suitable 
to apply a Monte-Carlo random walk, since the modified Green's functions used in each 
iteration are written in a separable form. We find with a slight modification that the 
iteration formula can be written in a form similar to the diffusion equation. There is, 
however, a crucial difference between them. In our new iteration formula, a higher order 
MCF wave function is generated from a lower order one. Therefore, the equation cannot 
be a simple iteration of the same function as in the diffusion equation. Second, in order to 
guarantee the orthogonalization of each MCF wave function, we must deal the difference of 
two MCF wave functions. This causes a sort of the "sign problem". [[□]] Here, we propose 



a correct method of counting the weight in the random walk procedure for an equation of 
a non- iterative type. As a result, we find out that our method is closely related to "the 
multiple cancellation" method. [0 With our cancellation method, we can take care of the 
sign problem in our formula. 

The validity of our method of counting weight and the convergence of the MCF are ex- 
tensively examined for a simple bound state problem in one-dimension including the excited 
states. As a more realistic example of a multi-dimensional problem, we study A = 4 nuclei 
with the Volkov potential. ]14|] Our method is successful in the treatment both examples. 

We briefly review the method of continued fractions in §2. Our random walk method 
in a non-iterative equation and multiple cancellation mechanism are discussed in §3. We 
present our results in §4, and a summary is given in §5. 
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II. METHOD OF CONTINUED FRACTIONS 



The method of continued fractions (MCF) has been developed as an efficient method 
to study scattering problems in atomic and nuclear physics. PJlQfl We briefly review this 
method, here. 

The Schrodinger equation of nuclear bound state can be written as 

|0) = G Q {E)V\<j>) , (1) 
where G (E) is a free Green's function, 

Go(E) = > ( 2 ) 

and Hq and V are the kinetic energy of the nucleons and the nuclear potential, respectively. 

In this paper, we use the modified Green's function (MCFG) approach, || since it is 
easy to implement in Monte-Carlo evaluations. We introduce the arbitrary functions \Fq) 
and \f) which are not orthogonal to the eigenfunction \<p). A new function \Fi) is defined as 

\F t ) = G (E)V\F ) , (3) 

and a modified Green's function G\(E) as 

G\(E) is chosen so as to operate only on the space which is orthogonal to V|Fo) as 

Gx(E)V\F ) = . (5) 
We eliminate Gq{E) from the Schrodinger equation ([!]) using (f|): 

\^=( G ±( E ) + Mnk) v \<t>) ■ (6) 



Then we define as 



and (|l|) can be written as 



(f\V\F ), 



Multiplying (§) by (f\V from the left, we obtain the following condition to determine the 
bound state energy E and the wave function |</>): 

Q(E) = (f\V\F ) - (f\V\<h) = . (9) 
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The bound state condition can be written in terms of a continued fraction as 

'^rw |fl> ' (10) 

|F m ) = G l (E)V\F l ) , (11) 

and 

G -< B » = G «< B > - jmh • < 12) 

where \<f>j) is the z-th order wave functions. The z-th Green's function Gi(E) satisfies the 
orthogonality relation 

Gi(E)V\Fj) = . (i=j + l) (13) 

We rewrite (|TUj) using (O), multiply (/|V from the left, and obtain the following iteration 
formula for \<pi): 

</M « = • (14) 

As a result, the bound state condition @ can be written in the form of a continued fraction: 

2 

=x — 2 . (15) 

x 2 

x i 72 — 



%2 - 



x 3 



■ T 3 - . 

Here, the matrix elements Xj are defined as 

Xi = (f\V\Fj . (16) 

In the MCF calculation, we only need to generate \Fj) and Gi(E) iteratively, following 
the iteration procedure (|ITD and (|T^). In this procedure, Gi(E) projects out the states 
which are orthogonal to the lower order wave functions, as indicated in (0). Therefore, Xi 
is expected to become smaller for larger i, and rapidly converges into zero. The advantage of 
this continued fraction series based on Ref. |L(| for the Monte-Carlo calculation is that the 



Xi are simply given by the matrix element of the potential between \Fi) and the analytically 
known function |/), in contrast to those given in Ref. 0, where we need the matrix elements 
between 

The bound state energy is determined by solving the equation 

Q(E) = (17) 

numerically. In actual practice, we perform a parameter search for E. Since we do not adopt 
any projections onto the ground state, we can also apply MCF to study the excited states 
without any restriction such as the e~ rH projection. 
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Once the energy of the bound state has been determined, the relevant eigenfunction is 
given as 

\(j))=Af[\F l )+y 1 \F 2 )+y 1 y 2 \F 3 )+y 1 y 2 y 3 \F 4 ) + ---} , (18) 
where is defined as 

Vi = x t - (f\v\M ' (19) 

and AT is a normalization constant. 



III. CONTINUED FRACTIONS AND MONTE-CARLO METHOD 

The original form of the MCF iteration formula flTT|) is not suitable for Monte-Carlo 
evaluations. The difficulty is caused by the Green's functions Gi(E). Since the Gi(E) obey 
the iteration formula (|T2|) and are of a separable form, 

Gl {E) = G^{E) , (20) 

•Ei—l 

it does not seem practical to generate the random numbers according to the probability of 
Gi(E). In order to resolve these problems, we use the orthogonality relation flT3|) and modify 
the iteration formula flTTp as follows: 



= G (E)V 



Xi-1 



(21) 



The advantage of this formula is that a higher order wave function \F i+ i) is always generated 
by the free Green's function Gq{E), which has a simple analytical form. Therefore, we can 
use the Monte-Carlo random walk technique in the MCF iteration procedure. 

However, it should be noted that the iteration formula ([H]) is not of the same form as 
the original Schrodinger equation. The Schrodinger equation, 

10) = Gq{E)V\4>) , (22) 



can be regarded as a diffusion equation. \(p) will converge to a ground state by iterating (|22|) 
with a Monte-Carlo random walk. In contrast, (|2lf) can be written as 

= G (E)V\^ ) , (23) 

where and \ipo) represent \F i+ i) and [\Fi) — \Fi^i)xi/xi-i], respectively. In our iteration 
formula (p3|) , the wave function has different functional form from |^o) an d is generated 
from the two terms \F{) and Therefore, we cannot naively apply the random walk 

technique of the diffusion equation. Otherwise, a destructive interference of and 



in (21) may cause the "sign problem". 

In this section, we extend the random walk technique for ( |2"B"D and propose a new method 
of counting weight. The destructive interference of two terms in ([H]) can be taken care of 
by our proper method of counting weight. 
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A. Random walk in MCF 



We study an algorithm to apply the random walk method with the free Green's function 
G (R, R!) to our following iteration formula of MCF: 

tp^R) = J G (R, R')W(R')MR')dR' ■ ( 24 ) 
Here we re-scale the variables so that the free Green's function Gq(R, R') is normalized as 



-V£ + 1 



G (R, R') = 5(R - R') , (25) 



where R denotes an M-dimensional vector, i.e. R = {z\, ■ ■ ■ , zm} defined in the standard 
manner given in Ref. [13]. The potential W(R) is defined as 



W(R) = -V(R)/B , (26) 

where B is the binding energy. 

For simplicity we assume that the wave functions ij)\(R) and i[)q(R) and the potential 
W(R) are positive definite. Note that the positive sign of W(R) implies that this is an 
attractive potential. We assume ip\(R) and W(R)tPq(R) are represented by the sum of 
points {Rk,i} distributed with the probability ipk(R) and the weights {Wfc,i}, where k = 0,1. 
The question is how {Rij} and {Wij} should be generated from given sets of {Roj} and 
{Wqj}. The relation can be obtained by studying the following overlap integral of ipi(R) 
with an arbitrary function F(R): 

1 = J F{R)%l)i{R)dR . (27) 

X is given by using {Rij} and {W^ij} as 

N 

l = T, W hi F (Ri,j) , (28) 
j 



while the same integral is expressed by {Roj} and {Woj} using (03) as 

N 

1 = J F(R) W ,iGo(R, Ro,i)dR , (29) 

i 

where N is the number of sampling points. 

We multiply the integrand of fl2"9"|) by unity in the form Gq(R, X)/Gq(R, X) and obtain 



1 = J A(R,X)G (R,X)dR , (30) 

where A(R, X) is defined as 

A(R, X) = F(R) £ W f^ R R ^ . (31) 
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Here, X can be an arbitrary vector in M- dimensional space. We impose hyper spherical 
boundary conditions on Gq(R,X) for the bound state problem, 



G (R,X) = for \R\,\X\^oo. (32) 

Gq(R,X) becomes a function of distance \1Z\, where 1Z is defined as 71 = R — X. Since 
Gq(R,X) is normalized as 

j G (R,X)dR = J G (\K\)dK = 1 , (33) 

Go(|72.|) can be regarded as a probability function. 

We change the integral variable from R to 1Z in (|30|) and generate a set of random 
numbers {TZj} with probability Gq(\7Z\). The algorithm to generate {TZj} in M dimensional 
space is given in Ref. fl3| . Equation ( j30|) is rewritten as follows: 



X -\ F{ X + Uj) NGoix + n^x) ■ (34) 

Since Q3"3| ) holds for any values of X, we choose X as Rqj for the j-th random number TZj. 



Comparing fl3"4|) with (|28|), we can define the new points Rij and the weight Hij as 

Ri,j = Ro,j + TZj , (35) 

and 

w ^ _ E» W / o^G (-Roj + 7^.j, R ,j) 
' NG (R j + TZj, Ro t j) 

As a result, fl3l>D and ( |3T)D show how the distribution of points {Rij} and the weights {W^ij} 
of ipi(R) should be generated from those of W (R)ipo(R) . A new point R\j is generated by 
a random walk from Rqj similar to the diffusion equation. However, a new weight Wij is 
generated by taking into account all of the old points {-Ro,i} and weights {Ho,?} through the 
free Green's function. In the naive random walk method, Wij is simply given by Ho,i, which 
projects into the ground state wave function. For a non-diffusion type integral formula like 
([24]), weights of new points should be generated by using this method of counting weight. 



B. Multiple cancellation 

In general, ipo(R) and W(R) are not positive definite, and we must include their sign 
information in our formula. This can be done by adding the sign sqj to the weight Wqj as 

W 0J > 8 jW 0J • (37) 

Similarly, adding the sign s^j to Wij, fl36|) is replaced by 

s _ Ej so ! iWo,iGo(R j + TZj, Rpj) 

' J ' NG (Rqj + TZj, Rqj) 
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Equation ( p8| ) is quite similar to the "multiple cancellation" method proposed by J. B. 
Anderson et al.. [0 The multiple cancellation method with a random walk has been applied 
to atomic physics, where it has provided very accurate results. In the conventional random 
walk technique, points with positive and negative signs walk independently. Then the weight 
of each point grows, and a destructive interference of these wave functions induces large 
numerical errors. On the other hand, our formula (j38|) shows that weight and sign of new 
points are generated by taking into account all the weights and signs of old points. Therefore, 
the points with positive and negative weights are canceled out in the course of counting the 
new weight. Equation (581) plays an essential role in the MCF calculation, and its importance 



will be demonstrated in the one-dimensional example given in the next section. 

IV. APPLICATIONS TO BOUND STATE PROBLEMS 

The procedure of our method is summarized as follows: 
1. Assume an appropriate form of the trial function \Fq). 



2. Distribute the points by the Metropolis method [T5j according to the functional form 
of|F ). 

3. Calculate the matrix element X{ of the i-th iteration step. 

4. Generate new points of by a random walk according to the procedure described 
in the previous section. 

5. Iterate steps 3 and 4 until the X{ and Q(E) converge. 

The binding energy B(= —E) is determined by searching for an energy that satisfies the 
condition Q(E) = 0. 



A. One-dimensional model 

As a simplest example, we study a one-dimensional potential model. In this one- 
dimensional problem, we can study the validity of our method by comparing with a conven- 
tional numerical method. We take the mass of the particle as m = 940 MeV and a simple 
central potential V(x) of exponential form with a one pion range as follows: 

V(x) = V exp(-a\x\) (39) 

with 

V = -69.3 (MeV) , a = 0.711 (far 1 ) . (40) 

The strength of the potential is chosen so as to produce three bound states at —38.7, —10.6 
and —2.01 MeV. The binding energy of the second excited state is small, and is compared to 
the deuteron binding energy. Since this model is one- dimensional and does not contain any 
other degrees of freedom, such as spin or isospin, the wave function of the bound states can 
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be classified by their parity. The ground and second excited states are even parity states, 
and the first excited state is an odd parity state. The validity of our method for the second 
excited state is particularly interesting, since it has the same symmetry as the ground state. 
Introducing a dimensionless variable z defined as 



z = V2mBx , (41) 

the free Green's function is given by 

G (z,z') = ±ex V (-\z-z'\). (42) 

The trial function is taken to be a Gaussian form as follows: 

1 x 2 

\Fo) = \f) = i — ex p(~^) for even P arit y> ( 43 ) 

J by/% 10 



I 2 x x 2 
= V b eXp( ~2&2 ) for odd parity - (44) 

We assume 6 = 3 (fm) for the ground and first excited states and b = 10 (fm) for the second 
excited state. 

In the Monte-Carlo calculation, the numbers of sampling points were 20,000, 70,000 and 
100,000 for the ground, first, and second excited states, respectively. We needed a larger 
number of points for the excited states than for the ground state in order to obtain reasonable 
accuracy. The initial distribution of points with the probability ||-F )| was generated by the 
Metropolis method. A typical number for the thermalization step was 100 in this problem. 

The fast convergence of the MCF iteration is related to the orthogonality relation in 



(fL3[), which is implemented in our iteration formula in (pT|). As an example, we consider the 



second order MCF wave function |F 2 ), 

X 1 



G (E)V 



•i'o 



(45) 



\F 2 ) is generated by two terms, \F{) and — |i<o)a;i/a;o, which will cancel each other. The higher 
order \Fi) decreases quickly with this cancellation mechanism. Therefore, accurate treatment 
of this cancellation is important. If we generate |_F 2 ) from two terms independently, the 
resultant distribution of the points for \F 2 ) is given by the large overlap of the points with 
positive and negative signs shown in Fig. [I] at the energy of the second excited state. \F 2 ) 
generated from these points will include large errors, which is just our "sign problem". 
Because of this sign problem, it is difficult to obtain convergence for the higher order MCF 
wave functions in the standard random walk method. On the other hand, using our method 
of counting weight in (^), the cancellation between two terms is performed analytically, 
and the distribution of the positive points for \F 2 ) is clearly separated from the negative 
points, as shown in Fig. 0. Our method of counting weight plays an important role in MCF 
iteration with the Monte-Carlo method. 

MCF wave functions at each iteration step are shown in Fig. |3|. They decrease rapidly 
as the iteration steps increase, even for the second excited state. This agrees well with the 
results obtained using the standard method of numerical integration. 
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Consequently, the matrix elements Xi defined in (|16D and Q(E) also converge rapidly as a 
function of the iteration step as shown in Fig. [| and |5], respectively. A few steps of iteration 
is enough to obtain reasonable convergence for the ground and first excited states. Slower 
convergence for the second excited state is due to our simple Gaussian trial function, even 
though one expects nodes in the wave function of the second excited state. 

The energy dependence of Q(E) is shown in Fig. |6| for the ground state. To determine 
the binding energy from the condition Q(E) = 0, we fitted Q(E) by using the linear energy- 
dependence of Q(E) = near the binding energy with the least square search method. Since 
values of Q(E) fluctuate due to sampling errors, we performed calculations several times with 
different sets of random numbers. Averaging these results, we obtained the energies of three 
bound states and errors estimated by the standard deviation, as shown in Table |. Our 
results with the Monte-Carlo method are in a good agreement with those of the 'exact' 
calculation for all three states. 

Once we obtained the binding energies, the eigenfunctions were given by the MCF wave 
functions \Fi) and matrix elements Xj. The wave functions for the ground, first, and second 
excited states are shown in Fig. |7|. The wave functions are also in good agreement with the 
'exact' calculation. Note that we started with a trial function without nodes for the second 
excited state. We could use more realistic trial functions with nodes, and if we did so we 
would obtain faster convergence. However, our result shows that it is not necessary to know 
the correct nodes of the wave function when we construct a trial function in this method. 
In summary, our new Monte-Carlo method works for the calculation of binding energies and 
wave functions including the excited states, and we believe that our method will be useful 
in studying few-nucleon bound state systems. 



B. Four-nucleon bound state as an example of a multi-dimensional problem 

As a realistic example, we study a four-nucleon bound state in three-dimensional space. 
The Schrodinger equation is expressed as follows: 

-E^-V^|0)+E^i0) = -^)- (46) 

k=i i<j 

Here, the are the Jacobi coordinates defined as 

' T\ = X 2 ~ X 1 , 

r 2 = x 3 - (x 1 + x 2 )/2 , 
r 3 = x 4 - (x 1 + x 2 + x 3 )/3 , 

and the /i^ are reduced masses corresponding to the coordinates r k : 

/ii = m/2 , 
fi 2 = 2m/3 , 
/i 3 = 3m/4 . 

Here, we adopt the Volkov potential, which is a central potential and has a soft core at 
short range: |14[ 
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Vij = Vx exp(-air 2 ) + V 2 exp(-a 2 r 2 ) , (47) 

with 

(Vt = -83.34 (MeV) , a x = 1.6CT 2 (fm~ 2 ) , 
\V 2 = 144.86 (MeV) , a 2 = 0.82~ 2 (fin" 2 ) , . 

Since the ground state energy of 4 He using this potential has been studied by various methods 
such as the Faddeev-Yakubovsky (FY) method taking into account only the S"-wave com- 
ponent, |l] the stochastic variational method (SVM) Jl6| and the hyperspherical harmonics 
(HH), |T7| it is a good example to examine our method in a four-body system. 



The free Green's function in dimensionless coordinates is expressed as 

G ° Wfl ' )= 2(^( i+ ! + # + 1) exp( - K) • (48) 

where 1Z = \R—R'\ is a relative distance in the nine- dimensional space. As the dimensionality 
increases, the free Green's function becomes more singular at 1Z — 0. It is rather difficult to 
treat the singularity of the free Green's function with the conventional integration method. 
Since the random walk method inverts the free Green's function analytically, we can obtain 
numerically stable results. A standard method to generate points with the probability Gq(TZ) 



is given in Ref. [[13 



The ground state of 4 He is expected to be spatially symmetric and anti-symmetric in 
spin-isospin space. Then we take a trial function |F ) given by the product of a Gaussian 
wave function \ip s ), a correlation function g(r) and the spin-isospin function |H), as 

\F )=AfY[g(r l ,)\ iPs )\E) , (49) 

i<j 



with 



|^> = nexp(-r 2 /26 2 ) , (50) 

i<j 

where bo is a size parameter. We used the form of correlation function g(r) 

9(nj) = [1 - a exp(-r 2 /6 9 )] 2 (51) 

with 

a = 0.11 , b g = 0.5 (fm 2 ) . 

The parameters of the correlation function are determined to simulate the wave function of 
a two-nucleon bound state in this potential. The size parameter bo of the Gaussian wave 
function in (|50| ) is simply determined by the variational method with one parameter. For 
the trial function |/), we assume the same functional form as and take a slightly larger 
value of the size parameter bo, 

b = 3.0 (fm) : for |F ) (52) 
= 4.5 (fm) : for |/) . (53) 
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Since the Volkov potential is a central potential with no exchange terms, the symmetry of 
the spatial and spin-isospin wave functions are preserved separately. 

The convergence of Xj is displayed in Fig. |8|, which shows that three-times iterations are 
enough in this example. Similarly, Q(E) also converges rapidly, as shown in Fig. H Here, 
we take N = 200, 000 sampling points. 

The bound state energy and its uncertainty are estimated in the same way as in the one- 
dimensional case. We assume a linear function for fitting the energy dependence of Q(E), 



as shown in Fig. [10]. As a result, we obtain E = —29.88 ± 0.59 (MeV). Our results for the 
binding energy of 4 He are in good agreement with the results of the other methods within 
uncertainties, as tabulated in Table It is shown that our method also works for higher 
dimensional examples. 



V. SUMMARY AND DISCUSSION 

The Monte-Carlo method is useful in the study of the bound state of the many-nucleon 
systems because of its simplicity and efficiency for multi-dimensional calculations. In this 
work, we propose a new method to study the bound state of a few-nucleon system without 
using the diffusive-type projection method. Therefore it is expected to work not only for 
the ground state but also for excited states. 

We used the method of continued fractions (MCF) as the main formalism of the calcu- 
lation. The advantages of this method are its ability to calculate the excited states and the 
rapid convergence of the continued fraction series. We modified the original formula of MCF 
into a suitable form for the Monte-Carlo random walk method. We generated MCF wave 
functions by the random walk method. Here, it was shown that the cancellation mechanism 
of lower order MCF wave functions is important, and we should use our method of counting 
weight in this calculation to avoid the "sign problem" . 

Our method was examined in application to a one-dimensional problem in detail. We 
demonstrated that our new method works even for the calculation of excited states. We also 
studied the bound state of a four-nucleon system and obtained results that compare well 
with those of other methods. These facts suggest that our method is a powerful method to 
study the bound states of few-nucleon systems. 

Application of this method to few-nucleon systems with realistic nuclear potentials in- 
cluding exchange terms is one of our future problems. Here, we have to deal with the coupled 
channel problem of wave functions with various types of exchange symmetry. It is impor- 
tant to retain the correct symmetry of the wave function in each iteration step. It is also 
necessary to take into account the correlation of nucleons in the trial function as much as 
possible. The wave function obtained by a fully variational calculation should be used as a 
more plausible trial function in our procedure. This will also accelerate the convergence of 
the continued fraction series. 
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FIG. 1. MCF wave functions at energy (2.01 MeV) of the 2nd-excited state. The solid line, 
dotted line, and dashed line correspond to |i*2), \F\), and —\Fq)x\/xq, respectively. 
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FIG. 2. MCF wave function \F2) generated by our method of weight counting. The solid and 
dotted lines correspond to the population of the points for positive and negative weight, respectively. 
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X (fin) 

FIG. 3. The wave function of MCF, at the energy of the 2nd-excited state (-2.01 MeV). 
The solid lines are obtained by Monte-Carlo method, and the dashed lines are obtained by a 
numerically 'exact' method. The numbers (0 ~ 5) specify the wave functions \Fi). Normalization 
of the wave functions obtained by the Monte-Carlo method is adjusted to the 'exact' ones. 
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FIG. 4. Convergence of function of the iteration step i for the ground, 1st-, and 

2nd-excited states in units of MeV. Cross points were obtained by the Monte-Carlo method, while 
the box points were obtained by a numerically 'exact' method. 
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FIG. 5. Convergence of Q(E) as a function of the iteration step i at the ground state energy 
with arbitrary scale. The points marked by crosses and boxes were obtained by the Monte-Carlo 
method and a numerically 'exact' methods, respectively. 



18 



X10 



•1 



x Monte-Carlo ^ 



Linear-fitting 



x 



n(E) 



38 38.2 38.4 38.6 38.8 39 39.2 39.4 



Binding energy : B=-E (MeV) 

FIG. 6. Energy dependence of Q(E). Points marked by crosses were obtained by the Monte-Carlo 
method. The dotted-dashed line was obtained by a least square fitting. The units of the vertical 
axis are the same as in Fig. 5. 
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FIG. 7. The eigenfunction for the ground, 1st-, and 2nd-excited states. The solid lines were 
obtained by the Monte-Carlo method. The dashed lines were obtained by a numerically 'exact' 
method. Normalization of the wave functions obtained by the Monte-Carlo method were adjusted 
to the exact ones. 
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FIG. 8. Convergence of xi as a function of the iteration step i at B = 30 (MeV). 



21 








Iteration steps : i (times) 



FIG. 9. Convergence of Q(E) as a function of the iteration step i at B =30 (MeV). 
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FIG. 10. Energy dependence of £l(E), fitted using the same method as in the one-dimensional 
model. 
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state 


'exact' 


Monte- Carlo 


ground 


-38.7 


-38.8 ± 0.1 


lst-excited 


-10.6 


-10.6 ± 0.1 


2nd-excited 


-2.01 


-2.03 ± 0.14 



TABLE I. The eigenvalues (MeV) for the ground, 1st-, and 2nd-excited states. The results 
obtained by the Monte-Carlo method are compared with those obtained by a numerically 'exact' 
method. 
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Method 


Energy (MeV) 


Faddeev-Yakubovsky [1 


-30.27 


HH [ig 


-30.40 


SVM H 


-30.42 


MCF with MC 


-29.88 ± 0.59 



TABLE II. The bound state energy of 4 He. 



25 



